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We propose a gravity dual of antiferromagnetic quantum phase transition (QPT) induced by 
magnetic field and study the critical behavior around the quantum critical point (QCP). It turns 
out that the boundary critical theory is a strong coupling theory with dynamic exponent z = 2 and 
that the hyperscaling law is violated and logarithmic corrections appear near the QCP. Some novel 
scaling relations are predicated, which can be tested by experiment data in future. We also make 
some comparison with experimental data on low-dimensional magnets BiCoPOs and pyrochlores 


Er2—22 Y 2s Ti2O7. 


Introduction. —Quantum phase transition (QPT) and 
the behavior of quantum systems in the vicinity of the 
corresponding quantum critical point (QCP) have at- 
tracted a lot of attention both in theory and experi- 
ment sides recently [1-3]. In contrast to their classical 
counterparts induced by thermal fluctuations arising at 
finite temperature T > 0, QPTs happen at zero tem- 
perature and are governed by quantum fluctuations as- 
sociated with the Heisenberg uncertainty and are driven 
by a certain control parameter rather than temperature, 
e.g., composition, magnetic field or pressure, etc. In con- 
densed matter physics, such a quantum criticality is con- 
sidered to play an important role in some interesting phe- 
nomena [4, 5]. 


One of intensively discussed QPTs is ordered- 
disordered QPT in antiferromagnetic materials induced 
by magnetic field (see for example, Refs. [6—10]), espe- 
cially in the heavy-fermion systems, since they can be 
tuned continuously from an antiferromagnetic (AF) state 
to a paramagnetic (PM) metallic state by varying a single 
parameter [4]. In these materials, QPT naturally belongs 
to the phenomenon involving strongly correlated many- 
body systems [11-13]. However, the complete theoreti- 
cal descriptions valid in all the energy (or temperature) 
region are still lacking. In order to study and character- 
ize strongly coupled quantum critical systems, some new 
methods are called for. 


Thanks to the feature of the weak/strong duality, the 
AdS/CFT correspondence provides a powerful approach 
to study such strongly coupled systems. This duality 
relates a weak coupling gravitational theory in a (d+ 1)- 
dimensional asymptotically anti-de Sitter (AdS) space- 
time to a d-dimensional strong coupling conformal field 
theory (CFT) in the AdS boundary [14-16]. In recent 
years, we have indeed witnessed that the duality has 
been extensively applied into condensed matter physics 
and some significant progresses have been made [17-20]. 
In Ref. [21] we realized the ferromagnetic/paramagnetic 
phase transition in a holographic setup, and in Ref. [22] 


the holographic antiferromagnetic/paramagnetic phase 
transition was studied. We showed that the antiferro- 
magnetic transition temperature Ty is indeed suppressed 
by an external magnetic field and tends to zero when the 
magnetic field reaches its critical value Be. In this way 
the antiferromagnetic QPT induced by magnetic field is 
realized. However, it was shown that the model proposed 
in Ref. [21] contains a vector ghost, very recently, a mod- 
ified model was proposed [23], which is shown not only 
ghost free, but also causal well-defined, while it keeps 
the main results in the original model qualitatively. Here 
we will elaborate in some detail this QPT and study the 
corresponding critical properties in this new model. 
Holographic model. —In order to describe the spontaneous 
staggered magnetization which breaks the time rever- 
sal symmetry, we introduce two real antisymmetric ten- 
sor fields coupled with U(1) Maxwell field strength [22]. 
Based on the discussions in Refs. [22, 23], we take the 
bulk action as follows, 
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Here * is the Hodge star dual operator and dM denotes 
the exterior derivative of M. L is the radius of AdS 
space, 2k? = 167G with G the Newtonian gravitational 
constant, k, m? and J are all model parameters with 
J <0, à? characterizes the back reaction of the two po- 
larization fields M, (a) to the background geometry, and 
Liz describes the interaction between two polarization 
fields. Note that by rescaling the polarization fields and 


the parameter J, A? can also be viewed as the coupling 
strength between the polarization fields and the back- 
ground Maxwell field. In the AdS/CFT duality, the 
model parameters m and k are related to the dual opera- 
tor dimension in the boundary field theory, J is related to 
the self-coupling coefficient of the magnetic moment and 
k describes the interaction between two kinds of magnetic 
moments. The reason to introduce two antisymmetric 
fields for describing the anti-ferromagnetism was elabo- 
rated in Ref. [22]. Note that the form of v(m) is not 
unique, we choose this form as it can lead to spontaneous 
symmetry breaking (see Fig.1 in Ref. [23]) and to simplify 
the equations of motion of the model. Compared with 
the original model for antiferromagnetism in Ref. [22], 
the key change is to replace the covariant derivative of 
the polarization field M by the exterior derivative. This 
change can avoid the problems such as ghost and causal 
violation, while keep the significant results in the origi- 
nal model qualitatively and in addition this model has a 
potential origin in string/M theory [23]. 

In the probe limit of A — 0, we can neglect the 
back reaction of the two polarization fields on the back- 
ground geometry. The background we will consider is a 
dyonic Reissner-Nordstr6m-AdS black brane solution of 
the Einstein-Maxwell theory with a negative cosmologi- 
cal constant, and the metric reads [24] 
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Here both the black brane horizon r, and AdS radius L 
have been set to be unity. The temperature of the black 
brane is 


T = (3— p? — B?) /4r. (4) 


For the solution (3), the corresponding gauge potential 
is A, = a(1 — 1/r)dt + Bady. Here p is the chemical 
potential and B can be viewed as the external magnetic 
field of the dual boundary field theory. 

We consider a self-consistent ansatz for the tensor fields 
with nonvanishing components Mè”, May (a=1,2) and 
define 
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By this definition, the antiferromagnetic order parame- 
ter, i.e., the staggered magnetization, can be expressed 
as [22, 23] 


N'/X =- / drB/r?. (6) 


Then the antiferromagnetism phase corresponds to the 
case when Nt Æ 0. In our model, it just corresponds to 
the case of 8 £0, while a = 0. 


With this ansatz, it is found that the equations for 
Me ) are algebraic ones and can be solved directly [25]. 
Therefore we pay main attention on a and 8. At the 
horizon, the regular initial conditions should be imposed. 
The behavior of the solutions of equations in the UV 
region (near the AdS boundary) depends on the value of 
m? +k. When m? + k = 0, the asymptotic solutions will 
have a logarithmic term, we will not consider this case 
here. Instead when m? + k 4 0, we have the asymptotic 
solution as [25] 
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where a+ and (+ are all finite constants. To make the 
system condense into the antiferromagnetic phase, as in 
Ref. [22], we require that the term associated with the 
magnetic field B in Eq. (7) should be the leading term. 
When B = 0, we require that the condensation for 8 
appears spontaneously. With those, the parameters have 
to satisfy m? > k > 0 and 


Ji (k, m?) < J < Je (k, m’), (8) 


with JE (k, m?) = —(m? + k)’ (m? + 3/2 + k)/12 and 
a, = 6, = 0 according to the AdS/CFT dictionary (for 
details please see [25]). 

QCP, energy gap and spectrum.—Let us first consider 
the influence of the external magnetic field B on the anti- 
ferromagnetic critical temperature Ty. Near the critical 
temperature, the staggered magnetization is very small, 
i.e., 8 is a small quantity. In that case we can neglect the 
nonlinear terms of 6 and obtain the equations for a and 
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Here mé,o¢ and m3 g are two functions of a [25]. Without 
loss of generality, we can set 3(r;,) = 1. With increasing 
the magnetic field B from zero, the effective mass square 
of 6 increases, so that the critical temperature Ty de- 
creases. The critical temperature is plotted as a function 
of the external magnetic field in Fig. 1. When Ty is 
decreased to zero, an AdS2 geometry emerges near the 
horizon [25]. The existence of a stable IR fixed point in 
the emergent AdS2 region demands 
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at the horizon r = rp = 1. Then we can see that in 
the case of T = 0, when |B| < |B-|, 6 is unstable near 
the horizon and the condensation happens so that the 
staggered magnetization is no longer vanishing. When 
|B| > |B.|, however, 8 is stable at the horizon and the 


staggered magnetization is zero. Therefore, a QPT oc- 
curs at |B| = |B,| and the system is quantum disorder 
when |B| > |B,|. 

In order to investigate the magnetic fluctuations in the 
vicinity of QCP, we need to consider the perturbations 
of two polarization fields. To make the system be self- 
consistent at the linear order, the perturbations for all 
components of the polarization fields have to be consid- 
ered, 
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Put this perturbations into the equation of motions and 
compute to the 1st order for e, we can get their equa- 
tions of the perturbations (the details for perturbational 
equations can be found in Ref. [25]). In general, be- 
cause of the nonlinear potential, all the components 
of two polarization field couple with each other. Let 
B= (co) -08 )/2. In the paramagnetic magnetic phase 
(T > Ty or B > Be), when q,w —> 0, the equations for 
B decouple from others [25]. By imposing the ingoing 
condition at the horizon, it has the following asymptotic 
solution in the UV region, 
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According to the dictionary of AdS/CFT, up to a positive 
constant, the retarded Green’s function for 8 reads 


Gop = B_/B+. (13) 


Using the retarded Green’s function, we can define spec- 
trum function as P(w, 7) = ImG(w, 7)/7. When we 
turn on a small momentum 7, the energy of long-life 
quasi-particle, which corresponds to the peak of P(w, 7d), 
can be given by following dispersion relation, 


ws =A +e, €= =0. (14) 


Here A is the energy gap of quasi-particle excitation. In 
the vicinity of QCP, for the case of w = 0, the retarded 
Green’s function usually has the form of G ~ 1/(q? + 
1/€7), where € is called correlation length. At QCP, in 
general, the energy gap vanishes. Thus we have w, = €7. 
In addition, for small frequency and wave vector, we can 
define the dynamic exponent z in the way as wą ~ q”. 

Numerical Results. As the equations involved here 
are nonlinear, we have to solve them numerically. The 
different parameters satisfying restrictions (8) give sim- 
ilar results, we here therefore just take parameters as 
m? = 1,k = 7/8 and J = —0.67 as a typical example in 
the left plot of Fig. 1, and Figs. 2 and 3 and as J = —0.71 
in the right plot of Fig. 1 

We can see from Fig. 1 that the Néel temperature Ty is 
suppressed by external magnetic field. There is a critical 
magnetic field for given parameters, at which Ty is zero 
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FIG. 1. The antiferromagnetic critical temperature Ty versus 
the external magnetic field B compared with experimental 
data. The critical temperature Ty is calculated with the aid of 
the solution of Eqs. (9). At each fixed value of magnetic field 
B the temperature is determined at which the AF condensate 
starts to appear. Left: Comparion with experimental data 
from BiCoPOs Right: Comparion with experimental data for 
pyrochlore compounds: Er2—22 Yer TizO7. The experimental 
data are from [26, 28] and rescaled. 


and QPT occurs. After rescaled, different parameters 
give similar behaviors with some slight differences. The 
physical picture of this QPT can be understood as fol- 
lows. When magnetic field reaches its critical value, there 
magnetic spins become partially aligned along the direc- 
tion of the magnetic field. Therefore the system requires 
less thermal energy to destroy the remaining magnetic 
spins order. 

The holographic model can give some interesting scal- 
ing relations near the QCP. For small B, numerical re- 
sults show that Ty — Tyo «x B?, where Tyo denotes the 
critical temperature in the case without external mag- 
netic field. When magnetic field is close to Be, we find 
that Néel temperature is fitted well by following relation 


Tyn/ln Ty x (1— B/B), (15) 


where Ty = Ty /Tno. We will analytically present the 
relation (15) by considering the emergent geometry Ad S2 
in the IR limit [25]. 

When magnetic field B is larger than the critical value 
Be, the antiferromagnetic phase disappears even at zero 
temperature. In this case, the system comes into quan- 
tum disordered phase at zero temperature, in which there 
is a gapped magnetic excitation. In the left plot of Fig. 2, 
we show ImG with respect to the frequency of antiferro- 
magnetic excitation in the case with different magnetic 
field. In the case of 0 < B/Be— 1 « 1, there is a 
distinct peak which gives the energy gap for the exci- 
tation. With increasing magnetic field, the peak moves 
towards higher energy and becomes more and more indis- 
tinct. This means that the gap increases but the lifetime 
decreases when magnetic field increases. At the critical 
magnetic field B = Be, we see w = 0, which corresponds 
to a gapless long-lifetime antiferromagnetic excitation. In 
the region of B/B, — 1 — 0*, we find the energy gap is 
fitted well by following equation (see the right plot of 
Fig. 2) 


A x (B/Be — 1), with A = A/Tyo. (16) 
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FIG. 2. Left: The antiferromagnetic spectrum function in the 
case with different magnetic field when B > Be. Right: The 
gap energy versus the external magnetic field when B/Be — 
1> 0*. 
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FIG. 3. Left: GT! as a function of q? when w = 0 in the case 
of |1 — B/B-| = 0.01. The solid line is the fitting curve by 
G7! «x q’+1/é?. Right: The correlation length € versus the 
magnetic field when |B/Be — 1| > 0. 


In the left plot of Fig. 3, we plot the inverse Green’s 
function G35 (q) in the case of w = 0 and |1 — B/B,| = 
0.01. We can see that it obeys the behavior of G~! ~ q?+ 
1/€? as we expected before. Thus the Green’s function 
can give the correlation length by fitting the curve of G~+ 
as a function of q?, which is shown in the right plot of 
Fig. 3. We see that the correlation length € as a function 
of the tuning parameter B obeys the following relation 


Ex |B/B.-1|-”, with v ~ 1/2. (17) 


As to the dynamical exponent z, in antiferromagnetic 
metals, z = 2 [5]. In our holographic model, the dy- 
namical exponent can be calculated by using similar nu- 
meric method. The results indeed show that z ~ 2. 
This numerical result can be confirmed by the emergent 
AdS geometry in the IR region [25] and agrees with 
the result from a different holographic model proposed 
in Ref. [29]. Furthermore, from the energy gap (16), we 
see that this energy gap satisfies the universal scaling re- 
lation A ~ |B — B-,|*”, which further indicates z = 2 in 
this model. 

Discussion. -The relation we found in this paper of 
the Néel temperature with respect to the magnetic field 
in the vicinity of QCP is quite non-trivial. Note that 
the relation (15) is not the usual power-law behavior or 
square-root form. But it is an expected result in the 2-D 
QPTs in strong coupling case [5]. This non-trivial coin- 
cidence strongly implies a connection between these two 
different theories. In the antiferromagnetic metal where 


magnetic ordered is dominated by itinerant electrons, dy- 
namical exponent z is 2 near the QCP. Since our dual 
boundary theory is a 2-D theory, the effective dimension 
is thus deg = d+ z = 4, which is just the upper critical 
dimension of the Hertz field theory [2, 30]. In this case, 
the hyperscaling is violated and logarithmic correction 
behavior appears. In fact, the d = z = 2 quantum crit- 
ical theory is in general not a weak coupling theory at 
any T > 0. Instead, a strongly coupled effective classical 
model emerges that can be used to determine the critical 
dynamics [31]. Our results show that it can be described 
well by AdS/CFT correspondence and this provides a 
new example of the applicability of the gravity/gauge 
duality in condensed matter theory. 

It is quite interesting to apply this holographic model 
to realistic materials. Since the holographic model is 
independent of the microscopic details of the materi- 
als and their interactions, it should be suitable for a 
class of materials. Two potential AF-QPT materials 
are BiCoPOs; with critical magnetic field Be ~ 15.38T 
(which is obtained by fitting a power-law relation [26]) 
and Er2-97 Y2rTigO7 with Be ~ 1.5T for « = 0 [27, 28]. 
In Fig. 1 we present the experimental data for these two 
materials and the holographic results, where we choose 
the model parameters so that they can give out the best 
fitting. The holographic model gives Be ~ 16.2T and 
1.45T, respectively. While the experiment data show 
that the energy gap for Er2—2sY2xTi207 when B > Be 
obeys the linear relationship (16) with slope 4.2 (see Fig.8 
in Ref. [28]), our holographic model gives the slope 5.0 
with the chosen model parameters. These two slopes 
are in the same order. Note that these two kinds of 
material have different microscopic structures and com- 
plex interactions, it is remarkable that the simple holo- 
graphic model can give a self-consistent description for 
those two materials. In addition, it is also worth men- 
tioning here that the doping at the magnetic site (x up 
to 0.085) has a very little influence on the critical behav- 
ior of Er2—2xY2sTi207. This indicates that an emergent 
universal behavior appears in these materials from very 
different microscopic details and could be described by 
the holographic model. 

As the critical behavior of a QPT induced by a mag- 
netic field, the three scaling relations (15), (16) and (17) 
near the critical point are our main results from the holo- 
graphic model. Besides the energy gap (16), our predi- 
cations on the scaling relations (15) and (17) can also 
be confirmed by experiments. Unfortunately at the mo- 
ment they cannot be checked by the existing experiments 
because the experiment data are absent when the Néel 
temperature is very close to zero. Of course whether 
the holographic model is suitable for these two materials 
needs more evidence. It is also very interesting to find 
more materials satisfying the conditions of this model 
and to check our predictions. We expect that this model 
can be confirmed or falsified experimentally soon. 
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